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ABSTRACT 

O 
(N 

^' The evolution of primordial collapsing clouds and formation of proto- 

Population III stars are investigated using three-dimensional ideal MHD sim- 
ulations. We calculated the collapse of magnetized primordial clouds from the 
prestellar stage until the epoch after the proto-Population III star formation, 
spatially resolving both parsec-scale clouds and sub-AU scale protostars. The 
O ' formation process of proto-population III star is characterized by the ratio of 

w ■ rotational to magnetic energy of the natal cloud. When the rotational energy is 

"-^ ■ larger than the magnetic energy, fragmentation occurs in the collapsing primor- 

I ' dial cloud before the proto-Population III star formation and binary or multiple 

i-{ • system appears. Instead, when the magnetic energy is larger than the rotational 

52 ■ energy, strong jet with > lOOkms"^ is driven by circumstellar disk around the 

proto-population III star without fragmentation. Thus, even in the early uni- 
verse, the magnetic field plays an important role in the star formation process. 
> 

^ ' Subject headings: binaries: general — cosmology: theory — early universe — ISM: 

jets and outflows 

vn 

o 
o 

1. Introduction 

S^ . Magnetic fields is a key ingredient in present-day star formation. For example, proto- 

j^ . stellar jets, which are ubiquitous in star-forming regions, are considered to be driven from 

protostars by the Lorentz force. Protostellar jets influence gas accretion onto protostars and 
disturb the ambient medium. In addition, the angular momentum of the cloud is removed 
by magnetic braking and protostellar jets. The removal of angular momentum makes proto- 
star formation possible in a parent cloud that has a much larger specific angular momentum 
than the protostar. So far, magnetic effects in primordial gas clouds have been ignored in 
many studies because magnetic fields in the early universe are supposed to be extremely 
weak. However, recent studies indicate magnetic fields of moderate strength can exist even 
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in the early uni verse. Cosmologica l fluctuations produced magnetic fields before the epoch of 



recombination ( llchiki et al. II2006I ). These fields were sufficiently large to seed the magnetic 



fields in galaxies. A generation mecha nism for magnetic fields at the epoch of reionization 



are also proposed (JLanger et al. 
amplified up to ~ 10" 



I2OO3I ). in which magnetic fields in intergalactic matter are 
"^^ G. These fields can therefore increase up to ~ 10~^ — 10~^ G in the 
first collapsed object having number density of n ~ 10^ cm~^. These fields may influence the 
evolution of primordial gas clouds and formation of Population III stars. 

Under spherical symmetry including hydrodynar nical radiative transfer, many authors 
have carefully investigated both the present-day (e.g.. lMasunaga fc Inutsukal I2OOOI ) and pri- 



mordial (e.g.. lOmukai fc Nishilll998[ ) star formation processes. A significant difference be- 



tween present-day and primordial star formation exists in the thermal evolution of the col- 
lapsing gas cloud because of differences in the abundance of dust grains and metals. In 
present-day star formation, the gas temperature in molecular clouds is ~ 10 K. These clouds 
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collapse isothermally for ric ^ 10 cm~ . Then the gas becomes adiabatic at ric 
and an adiabatic core (or the first core) forms. After the dissociation of molecular hydrogen 
lO^^cm"^), the protostar forms at n^ — 10^^ cm~^. On the ot her hand, prirnordial 



[Ur 



gas clouds have ten iperatures of ~ 200 — 300 K at n^ ~ 10'^ cm ^ ( lOmukai et al. 1 12005 



Bromm et al. II2002I ). These clouds collapse keeping polytropic index 7 ~ 1.1 for a long 



range of 10^ cm ^ ^ ^c ^ 10^^ cm ^. Thus, the first core does not appear in the pri- 



mordial collapsing cloud. After the central density reaches ric — 10 cm , the thermal 
evoluti on of the primordial collapsing cloud begins to coincide with that of a present-day 



cloud ( lOmukai et al. II2005I ). The difference in thermal evolution between present-day and 



primordial clouds arises when ric ^ 10 cm . 

Another major difference between present-day and primordial star formation exists in 
their magnetic evolution. In present-day star formation, neutral gas is well-coupled with 
ions for ric ^ lO^^cm"^ and ric ^ lO^^cm"^, while the r nagnetic field dissipat es by Ohmic 



dissipation in the range of 10 cm ^ ^c ^ 10 cm ( iNakano et al. II2002I ). In the col- 



l apsing cloud, ~ 99% of the magnetic field is dissipated for 10 cm ^ ric ^ 10 cm 



( iMachida et al. 1120071 ). On the other hand, in a primordial gas cloud, the magnetic field 



couples strongly with the primordial gas during all phases of the star f ormation, as long as 
the i nitial field strength is weaker than Bq ^ 10"^(n/10^cm"^)°-^^ G (jMaki fc Susal |2004J . 
20071 ). In summary, the magnetic field is largely dissipated by Ohmic dissipation before pro- 
tostar formation in present-day clouds, while the magnetic field can continue to be amplified 
without dissipation in primordial clouds. 

The magnetic field in the collapsing cloud is closely related to the fragmentation or 
formation of binary or multiple stellar systems. In present-day star formation, the magnetic 
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field strongly suppresses rotation-driven fragmentation (JMachida et al. Il2008bl ). For primor- 
dial clouds, magnetic effects on fragmentation are still unknown. In this study, we investigate 
the evolution of weakly magnetized primordial clouds and the formation of Population III 
stars using three-dimensional simulations. 



Model and Numerical Method 



To calculate the evolution of primordial star-forming clouds in a large dynan iic range, we 
use the three-dimensional ideal MHD nested grid method (JMachida et al. ll2005al Jbl. l2006al Jbl). 
For gas pressure, we use a barotropic relation that approximates the result of one-zone 
calculation where the rmal and chemical processes for primordial gas are solved in detail 



flOmukai et al. II20051 ). 



As the initial state, we take a spherical cloud whose density is twice higher than that for 
hydrostatic equilibrium with external pressure (Bonnor-Ebert sphere). The central (number) 
density of the Bonnor-Ebert sphere is set at Pbe(O) = 3.5 x 10~^^gcm~^. The initial tem- 
perature is set at 250 K, and the radius of a Bonnor-Ebert sphere is Ro = 6.5 pc. To promote 
fragmentation, a small m = 2-mode density perturbation is imposed on the spherical cloud. 
The total mass contained inside Rq is M^ = 1.8 x 10^ Mq. The initial cloud rotates around 
the z-axis at a uniform angular velocity of Qq, and has a uniform magnetic field Bq parallel to 
the 2;-axis (or rotation axis). The initial model is characterized by two non-dimensional pa- 
rameters: the ratios of the rotational energy to the gravitational energy /3o (= -Erot/I-Egravl), 
and of the magnetic energy to the gravitational energy 70 (= -Emag/l-^gravl), where -Erot, -^mag 
and ^'grav are the rotational, magnetic and gravitational energies, respectively. We made 36 
models using different values for these two parameters. 



We adopted the nested grid method (for details, see iMachida et al. Il2005al . l2006al ) to 
obtain high spatial resolution near the center. Each level of a rectangular grid has the same 
number of cells (= 128 x 128 x 64), with the cell width h{l) depending on the grid level 
/. The cell width is halved with every increment of the grid level. The highest level of the 
grid changes dynamically: a new finer grid is generated whenever the minimum local Jeans 
length Aj falls below 8/i(/max), where h is the cell width. The maximum level of grids is 
restricted to /max = 30. Since the density is highest in the finest grid, generation of a new 
grid ensures the Jeans condition with a safety factor of 2. We begin our calculations with 
three grid levels (/ = 1 — 3). The box size of the initial finest grid / = 3 is chosen to be 2Ro, 
where Rq is the radius of the critical Bonnor-Ebert sphere. The coarsest grid {I = 1) then 
has a box size of 2^ Rq. 
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Results 



Figure [T] shows the cloud evolution before protostar formation [uc < 10^^ cni"'^) from 
the initial stage for magnetically-dominated model that has a larger magnetic energy than 
the rotation energy (70 > /3o). This model has parameters of (70, /3o) = (2 x 10~^, 10"''), 
and the initial cloud has a magnetic field strength of Bq = 10~^ G, and an angular velocity 
oi Qq = 2.3 X 10^^^ s^^. Figure [T]a shows the initial spherical cloud threaded by a uniform 
magnetic field. Figure [T]&-(i shows the cloud structure around the center of cloud when 
the central density reaches Uc = (b) 4.1 x 10'', (c) 7.0 x 10^°, and (d) 5.5 x lO^^cm"^, 
respectively. The density contours projected on the sidewall in these figures indicate that 
the central region becomes oblate as the cloud collapses because of the magnetic field and 
rotation, both of which are amplified as the cloud collapses. Figure [T]a-c? also shows that the 
magnetic field lines gradually converge toward the center as the central density increases. 
Figure Wg^i shows the evolution after the protostar formation, in which the strong jet is 
driven by the circumstellar disk. The speed of jet is ~ lOOkms"^ that corresponds to the 
Kepler speed near the protostar. 

Figure [2] shows a configuration of magnetic field lines (black-white streamlines) , a struc- 
ture of a jet (transparent iso-velocity surface) and fragments for some models against /3o 
- 7o plane. The magnetically-dominated models (70 > /3o) show a powerful jet after the 
proto-Population III star formation without fragmentation (or binary formation), while the 
rotation-dominated models (/3o > 70) show fragmentation and binary formation without jet 
driving. In the jet-driving models, the magnetic field lines are strongly twisted inside the jet 
region. These jets have hourglass-like configurations of the magnetic field lines, where the 
poloidal component is more dominant than the toroidal component. These configurations of 
the magnetic field lines can easily drive a strong jet by the disk wind mechanism. 

In fragmentation models, the configurations of the magnetic field lines are disturbed, 
and the toroidal field tends to be more dominant than the poloidal filed. Such configuration 
of the magnetic field lines means the weak magnetic field, which is insufficient to drive a 
jet. The weak poloidal field is attributed to the oblique shocks on the surface of the disk 
envelope, and the disturbance of the field lines are attributed to the orbital motion of the 
fragments and the spin of the fragments. The spin of the fragment winds the magnetic fields 
around its rotation axis, and the toroidal component of the magnetic field is amplified. Thus, 
in the further stages, the magnetic pressure may drive a jet. 
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4. Summary 

In this study, we calculated cloud evolution from the stage of ric = 10'^ cm~'^ until the 
protostar is formed (~ 10^^ cm~^) for 36 models, parameterizing the initial magnetic field 
strength and rotation, to investigate effects of magnetic fields in collapsing primordial clouds. 
Our calculations showed that fragmentation occurs but no jet appears when (3q > 70, and 
jet appears after the protostar formation without fragmentation when /3o < 70. Thus, in 
the collapsing primordial cloud, the cloud evolution is mainly controlled by the centrifugal 
force than the Lorentz force when /3o > 70, while the Lorenz force is more dominant than 
the centrifugal force when 70 > /3o- 

A jet is driven when the initial cloud has magnetic field of 

ric \2/3 



«»^i«nTnJ7^ (" W 



cm 



if the cloud rotates slowly as fi ^ 4 x 10^^^(nc/10'^ cm~^)^/^s~^. The power of a jet, e.g., 
a mass ejection rate, is considered to be controlled by the accretion rate as indicated in 
present-day star formation; the mass ejection rate of a jet is 1/10 of the mass accretion onto 
the protostar. The accretion rate of primordial star formation is expected to be considerably 
larger than that at present day, and it produces a stronger jet. The life time of the jet 
also seems to be controlled by accretion in the present-day; a jet stops when mass accre- 
tion stops. For Population I II stars, the gas accretion does not halt within their lifetimes 



(jOmukai fc Pallal I2OOII . l2003l ). Therefore, a jet also may continue during the all lifetime of 
the protostar, and the strong jet propagates to disturb a surrounding medium significantly. 
The disturbance of the medium could trigger the subsequent star formation as frequently 
observed in present-day star formation. 

Assuming the power law growth of Bq oc ric , the critical strength of the magnetic field, 
Bq = 10~^ G at ric = 10^ cm~^ corresponds to i?o = 5 x 10~^^ G at ric = 0. 01 cm~^ [see, 



eg. ([T|) ], which is much stronger than the background magnetic field derived by llchiki et al. 



( 20061 ). However, when the magnetic field is amplified to -B ~ 10 ^(ric/lO^cm ^y/^ G by 



some mechanisms (e.g.. ISchleicher et al.ll2010l ). the magnetic field can affect the collapse of 



the primordial cloud. Even if a cloud has a magnetic field weaker than the critical strength 
Bq = 10^^ G, the magnetic field may play an important role after the proto star formation. In 



analytical study of the evolution of accretion disks around the first stars (JTan fc Blackman 



20041 ) , it is suggested that magnetic fields amplified in the circumstellar disk eventually give 



rise to protostellar jets during the protostellar accretion phase. 

Rotation promotes fragmentation when the first collapsed objects has the angular ve- 
locity of fig ^ 10~^^(nc/10^cm~^)^/'^ s~^. The fragmentation is expected to produces binary 
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or multiple stellar system. When a multiple stellar system is formed, some stars can be 
ejected by close encounters. At the protostar formation epoch, the protostar has a mass 
of M ~ 10~^Mq. The ejected proto- Population III stars may evolve to metal-free brown 
dwarfs or low-mass stars. When a binary component in a multiple stellar system is ejected 
from the parent cloud by protostellar interaction, a low-mass metal free binary may also 
appea r in the early universe . It is considered t hat the extremely metal-poor ([Fe/H]<-5) 
stars (jChristlieb et al. II2OOII : iFrebel et al. II2005I ) are formed as binary members from metal- 
free ^as, and then have been polluted by the companion stars during the stellar evolution 
(ISuda et al. 1120041 ). In addition, a binar y frequency in Population III star niay be compara- 
ble to or larg er than that at present day (JKomiya et al. II2006I : lMachidall2008l : iMachida et al. 
2008d . l2009a| ]bl). In order to confirm the ejection scenario, the further long-term calculations 
are required. 
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Fig. 1. — (a-i) Time sequence of magnetically dominated model [{/3o, 70) = (10^'', 2 x 10^'^)] 
before protostar formation n^ < 10^^ cm"^. In each panel, the structures of the high-density 
region (isosurface) and magnetic field lines {black and white streamlines) are plotted in three 
dimensions, while the density contours {color and contour lines) and velocity vectors {thin 
arrows) are projected on each wall surface. The central number density ric, elapsed time t, 
grid level /, and grid size are also shown in each panel. The transparent red surface in panels 
{g)-{i) means the border between inflow and outflow inside which the mass ejected by the 
protostellar jet. 
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Fig. 2. — Final states in three dimensions against parameters /3o and 79. The magnetic field 
lines {black and white streamlines), high-density regions (isosurface), and jet [transparent 
isosurface) are plotted in each panel. The density contours [color and contour lines) and 
velocity vectors [thin arrows) are also projected on each wall surface. The grid level (/), and 
grid scale are shown in each panel. 



